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REMARKS 

The present invention is a method of active seismic monitoring of an 
underground formation as recited in claim 10 providing separation of induced 
microseismicity signals from seismic signals emitted during active seismic monitoring 
of an underground zone development and a method of active seismic monitoring of 
an underground formation as recited in claim 17 comprising discrimination of 
induced microseismicity signals from among signals emitted during active seismic 
monitoring of an underground zone under development. A method in accordance 
with the invention comprises carrying out seismic recording cycles with emission of 
seismic waves in the formation by coupling therewith at least one seismic source 
which, emits, simultaneously orthogonal signals so as to form a composite 
vibrational signal, receiving signals reflected by the formation in response to the 
emission of seismic waves, recording the signals received by at least one seismic 
pickup and processing the recorded signals to separate respective contributions of 
each seismic source to the received signals and to reconstruct seismograms 
equivalent to seismograms that would be obtained by separately actuating each 
seismic source, separating the induced seismic signals in the records from seismic 
signals resulting from active monitoring operations, by isolating a contribution thereof 
by comparison with a reference spectral model accounting for spectral contributions 
of each seismic source at emitted fundamental frequencies and at respective 
harmonics thereof, and by reconstructing the microseismicity signals by inversion in 
the time domain. 
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Claims 1-6 stand rejected under 35 U.S.C. §112, first paragraph, as not being 
enabled and further, claims 1-6 and 10 stand rejected under the second paragraph 
of 35 U.S.C. §1 12 as being indefinite pertaining to the recitation of a single source 
and multiple sources in the independent claims. 

Newly submitted claims 1 1-20, which correspond to claims 1-10, have been 
drafted to refer to "at least one seismic source" and thereafter refer to "each seismic 
source" which eliminates the issues of indefiniteness and insufficient disclosure. 
Additionally, the claims have been drafted to overcome the remaining stated 
grounds of rejection under the second paragraph of 35 U.S.C. §112 and have been 
further amended to improve their form for reexamination. 

Claims 1-10 stand rejected under 35 U.S.C. §101 as covering non-statutory 
subject matter "because the independent claims do not contain a tangible result, and 
therefore, do not pass the useful, concrete and tangible" test for statutory subject 
matter." The Examiner summarily concludes that the claims are drawn to 
mathematical operations involving the seismic and microseismic data and further 
concludes that the operations could be performed by a person in their head and 
therefore, they are not a tangible, real world result. These grounds of rejection are 
traversed for the following reasons. 

In the first place, the Examiner's contention that the claims recite mental acts 
is not supported by a review of the claimed subject matter or the specification. 
Specifically, claim 1 1 recites in the preamble a method of active seismic monitoring 
of an underground formation providing separation of induced microseismicitv signals 
from seismic sionals emitted within a context of active seismic monitoring of an 
underground zone and claim 1 7 recites In the preamble a method of active seismic 

8 



monitoring of an underground formation comprising discrimination of induced 
microseismicitv signals from among signals emitted within a context of active seismic 
monitoring of an underground zone under development. These claims recite actual 
seismic monitoring of an underground formation to provide separation of or 
discrimination of induced microseismicity signals which is the antithesis of the 
mental process and relates to processing of real world seismic signals to eliminate 
the effect of microseismicity signals. 

Claim 1 1 further recites at least one seismic source which emits 
simultaneously orthogonal signals so as to form a composite vibrational signal, 
receiving signals reflected by the formation in response to the emission of seismic 
waves, recording the signals received by at least one seismic pickup and processing 
the recorded signals to separate respective contributions of each seismic source to 
the received signals and to reconstruct seismograms equivalent to seismograms that 
would be obtained by separately actuating each seismic source which pertains to 
real world physical processes which cannot be performed mentally and have utility in 
monitoring underground formations, and claim 17 further recites carrying out seismic 
recording cycles with emission of seismic waves in a formation by coupling therewith 
seismic sources which emit simultaneously and are controlled by orthogonal signals 
so as to form a composite vibrational signal, receiving signals reflected by the 
formation in response to emission of seismic waves, recording of signals received by 
a seismic reception means and processing the recorded signals so as to separate 
respective contributions of the seismic sources to the received signals and to 
reconstruct seismograms equivalent to seismograms that would be obtained by 
separately actuating the seismic sources which, similar to claim 1 1 , pertain to real 
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world process operations that cannot be performed mentally and have utility in 
monitoring an underground formations. 

Finally, claim 1 1 further recites, in part, separating the induced 
microseismicity signals in the records from the seismic signals resulting from active 
monitoring operations, by isolating a contribution thereof by comparison with a 
reference spectral model by accounting for spectral contributions of each seismic 
source at emitted fundamental frequencies and at respective harmonics thereof, and 
by reconstructing the microseismicity signals by inversion in the time domain which 
are not mental steps and have utility in monitoring underground formations and 
claim 17 further, in part, recites calculating for each recording of a recording cycle 
respective contributions of the seismic sources at the fundamental frequency; 
calculating a ratio of a contribution to a current spectral model formed by updating 
the previous spectral model from frequencies emitted during a previous recording 
and from harmonics thereof; deducing a part of the recording n of a cycle p that can 
be associated with active seismic monitoring operations; deducing a part of the 
recording end of cycle p related to passive microseismic activities; forming by 
inversion seismograms that can be associated with active seismic monitoring 
operations by inversion in a time domain of the respective spectral contributions of 
each source at fundamental frequencies and at harmonics thereof, after completion 
of a measuring cycle; and forming underiying microseismic signals contained in the 
records by inversion in a time domain from a part related to passive microseismic 
activity which are not mental steps and have utility in monitoring underground 
formations. 
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Moreover, it is important to understand tliat the claimed separation or 
discrimination of induced microseismicity signals are signals resultant from active 
multisource seismic signals providing which are limited to monitoring operations in 
an underground zone under development, all of which are tangible, real world 
applications. See paragraphs 13-15 of the Second Substitute Specification for a 
discussion of the tangible application of the invention. 

Moreover, M.P.E.P. §21 06(b) defines a tangible result as "not necessarily 
meaning that a claim must either be tied to a particular machine or apparatus or 
must operate to change articles or materials to a different state or thing." In fact, the 
of separation of or discrimination of induced microseismicity signals from seismic 
signals emitted within a context of active seismic monitoring of an underground zone 
under development, as recited in claims 1 1 and 17, is a tangible result as defined 
under the foregoing section of the M.P.EP, 

The utility requirement defined under M.P.E.P. §21 06(a) is clearly met by 
claims 1 1 and 17 since it is understood by persons of ordinary skill in the art that 
microseismicity signals are in essence noise which are claimed as being separated 
out or discriminated from recorded seismic signals from multiple sources. It is 
understood that the best possible seismic information obtainable during monitoring 
operations in an underground zone under development as defined in paragraphs 
[0001] and [0013] of the Second Substitute Specification does not contain 
microseismicity signals. 

Finally, a concrete result under M.P.E.P. §2106 (c) is required to be 
reproducible. There is no doubt that the claimed processes of the claims recite the 
production of seismic signals which are free from microseismicity signals, which is 



the most accurate form of seismic signals, which facilitate active seismic monitoring 
operations. These signals meet the reproducible requirement of M.P.E.P. §2106(k). 

Claims 1-10 stand rejected under 35 U.S.C. §103 as being unpatentable over 
United States Patent 6,442.489 (Gendelman) in view of United States Patent 
5,721 ,71 0 (Sallas). These grounds of rejection are traversed for the following 
reasons. Even if the proposed combination of Gendelman in view of Sallas were 
made, the subject matter of the claimed invention as defined in independent 
claims 1 1 and 17 would not be achieved. 

Gendelman et al disclose a method for detection and monitoring of 
hydrocarbons which utilizes at least one receiver of seismic vibrations, a passive 
information signal is recorded in one embodiment, and in another embodiment, 
seismic vibrations are generated and an active information signal is recorded before 
and during the generation of the seismic vibrations. Neither embodiment separates 
or discriminates induced microseismicity signals from seismic signals resulting from 
active monitoring operations by isolating a contribution thereof of microseismicity 
signals by accounting for spectral contributions of each seismic source at emitted 
fundamental frequencies and at respective harmonics thereof and by reconstructing 
the microseismicity signals by inversion of the time domain as recited in claim 1 1 
and the sequence of steps a) - f) as recited in claim 17 which causes discrimination 
of induced microseismitic signals from among signals emitted within a context of 
active seismic monitoring of an underground zone under development. 

Gendelman's process uses two steps which first records microseismicity 
signals from the earth occurring prior to the activation of the seismic source and then 
emits an active signal in the earth and records the seismic signals which are 
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superimposed with the previously recorded seismicity signals. The microseismicity 
signals are detected by comparing frequency spectra of the signals received from 
the earth and then during application of seismic signals, resulting from active 
monitoring operations, isolates a characteristic spectral anomaly permitting the 
location of hydrocarbons. 

The present invention is substantially different from Gendelman in that there 
is no separation or discrimination of induced microseismicity signals from seismic 
signals emitted within a context of active seismic monitoring of an underground zone 
under development. 

Sallas et al do not cure the foregoing deficiencies. Sallas has been cited by 
the Examiner as teaching processing of data to reconstruct seismograms equivalent 
to those that would be obtained by actuating seismic sources separately as stated at 
the top of page 6 of the Office Action regarding the discussion of claim 1 and further, 
has been cited at the top of page 9 regarding the discussion of claim 7 that "Sallas 
teaches using multiple vibrators by hydrocarbon monitoring that emit orthogonal 
signals so as to form a composite vibrational signal... Sallas teaches processing the 
data to reconstruct the seismograms equivalent to those that would be obtained by 
actuating the seismic sources separately." However, if Sallas were combined with 
Gendelman utilizing the foregoing teachings of Sallas as relied upon by the 
Examiner, the foregoing separation of induced microseismicity signals from signals 
emitted within the context of active seismic monitoring of an underground zone 
under development including the claimed process steps, as recited in claim 1 1 , and 
discrimination of induced microseismicity signals from among signals emitted within 
a context of active seismic monitoring of an underground zone under development. 
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including the claimed process steps, as recited in claim 17, would not occur. The 
subject matter of claims 1 1 and 17 is not obvious since Sallas does not cure the 
deficiencies of Gendelman as discussed above. 

Moreover, dependent claims 1 2-1 6 and 1 8-20 define further aspects of the 
present invention which are not rendered obvious by the proposed combination of 
Gendelman and Sallas. 

In view of the foregoing amendments and remarks, it is submitted that each of 
the claims in the application is in condition for allowance. Accordingly, early 
allowance thereof is respectfully requested. 

To the extent necessary, Applicants petition for an extension of time under 
37 C.F.R. §1 .136. Please charge any shortage in fees due in connection with the 
filing of this paper, including extension of time fees, to Deposit Account No. 01-2135 
(612.44921X00) and please credit any excess fees to such Deposit Account. 



Respectfully submitted. 



ANTONELLI, TERRY. STOUT & KRAUS, LLP 




Donald E. Stout 
Registration No. 26,422 
(703)312-6600 
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SECOND SUBSTITUTE SPECIFICATION 
PATENT 

METHOD OF SEPARATING MICROSEISMICITY SIGNALS 
FROM SEISMIC SIGNALS EMITTED BY ONE OR MORE SOURCES 

Inventors : Julien MEUNIER and Frederic HUGUET 

ABSTRACT 

This invention is a method of separating induced microseismicity signals from seismic 
signals acquired within active seismic monitoring operations carried out in underground 
zones under development which has application for monitoring of underground 
hydrocarbon or fluid storage reservoirs. Seismic records are formed from signals emitted 
by one or more seismic sources controlled by orthogonal signals. In this case, the signals 
are processed to separate the respective contributions of the at least one seismic source to 
the signals received and to reconstruct the seismograms equivalent to those that would be 
obtained by actuating the at least one seismic source separately. The induced 
microseismicity signals (passive seismic signals) are separated from the seismic signals 
coming from sources (active seismic signals) essentially by isolating the contribution 
thereof by comparison with a spectral reference model at the fiindamental frequencies 
emitted and at the respective harmonics thereof, and by reconstructing in the time domain 
the induced microseismicity signals and the signals coming from the at least one seismic 
source. 



BACKGROUND OF THE INVENTION 
Field of the Invention 

[0001] The invention relates to a method for discriminating microseismic signals 

induced by the development of an underground reservoir, from among seismic signals 
emitted within the seismic monitoring operations. 

Description of the Prior Art 

[0002] In order to follow the evolution of an underground reservoir under 

development, it is well-known to permanently install seismic sources and seismic 
receivers, at the surface or in wells, whether for reservoir exploration or for development 
wells, and to carry out, at predetermined time intervals, recording of the seismic signals 
from the formation. These signals may be seismic signals induced by the site development 
activity or reflected by the geologic discontinuities of the formation in response to the 
emission of disturbances in the ground by means of one or more seismic sources. 
Comparison of the recordings performed at different times of the site development 
provides information on the evolution thereof. 

[0003] Various long-term seismic monitoring systems are described for example in 

the following patents EP-591,037 corresponding to US patent 5,461,594, French Patent 
2,593,292 corresponding to US Patent 4,775,009, and French Patent 2,728,973 
corresponding to US Patent 5,724,3 II) or French Patent 2,775,349, 

[0004] EP Patent 748,457 corresponding to US Patent 5,724,311 filed by two of 

the applicants describes a permanent monitoring system allowing obtaining perfect 
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reproducibility of the operating conditions in repetitive seismic monitoring operations 
carried out in an underground zone traversed by at least one well or borehole, notably in 
an underground gas storage reservoir. This system comprises, installed on a permanent 
basis, one or more set(s) of seismic receivers (buried at the surface or placed in one or 
more wells), several repetitive seismic sources (buried or at the surface), and a permanent 
linking network for selective energy supply to these sources. A central station remotely 
controls selectively each one of these sources and records the seismic signals coming from 
the underground zone in response to the seismic waves transmitted selectively in the 
ground by the sources. 

[0005] All these permanently installed sources whose coupling with the 

surrounding formations remains stable and this supply network, at least partly buried and 
whose surface coverage area is limited, allows carrying out a series of long-term seismic 
monitoring operations under stable operating conditions, without any risk of 
incompatibility with the development site activities. 

[0006] French Patent 2,775,349 corresponding to US Patent 6,182,012 filed by the 

same two applicants describes a method and a device intended for permanent monitoring 
of an underground formation in which cavities are formed for seismic sources (preferably 
deep enough to be coupled with the formation below the weathered zone (WZ)), at least 
one well is drilled either below each source or in the immediate vicinity thereof An 
antenna having a plurality of seismic receivers (geophones, hydrophone, etc.) is lowered 
into each well thus drilled and coupled with the surrounding formation. With this device, 
emission-reception cycles are carried out with emission of seismic waves in the formation 
by at least one source and acquisition of the signals picked up by the seismic receivers, in 
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response to the waves reflected by the formation. The receivers can for example be 
associated with a fluid transfer tube lowered in the well in order to connect an 
underground zone to a surface development device and coupled with the formation 
surrounding the well. 

[0007] The various sources of the seismic monitoring device can be actuated 

successively with a sufficient time interval between the triggering times for reception of 
the waves reflected by the zone being investigated. It is also possible to use seismic 
sources emitting the same signals that are triggered simultaneously to increase the power 
emitted. 

[0008] French Patent 2,805,051 filed by the applicants describes a method and a 

system for seismic monitoring of an underground formation. As diagrammatically shown 
in Figs.l to 4, the system comprises for example a network of seismic antennas 2 each 
having a series of seismic pickups 4 arranged at regular intervals along a well 3 drilled in 
the ground. This network can be regular as shown in Fig.2 or irregular. The pickups are, 
for example, vertically oriented one-directional or multi-axis (triphones) geophones and/or 
hydrophones. A seismic source 5 is arranged near to each antenna 2. Piezoelectric type 
vibrators such as those described in French Patent 2,791,780 US Patent 6,338,394 filed in 
the joint names of the applicants are advantageously used as sources and permanently 
installed in the immediate vicinity of each antenna 2. An electronic control and recording 
module 6 (Fig.l) can be associated with each antenna. These electronic modules 6 can be 
connected to a central control and synchronization laboratory 8. According to another 
implementation mode diagrammatically shown in Fig.4, antennas 2 are directly connected 
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to a central laboratory 8 including all the individual electronic control and recording 
modules 6 (Fig.l). 

[0009] Seismic waves are emitted in the underground formation by one or more 

seismic sources (vibrators) 5. In the case of multiple sources, the latter emit 
simultaneously and are controlled by orthogonal signals so as to form a composite 
vibrational signal. The vibrators work permanently or at least for a sufficiently long time 
to obtain a sufficient signal-to-noise ratio. The seismic waves generated by seismic 
sources 5 are propagated downwards (downgoing waves 9). These incident waves are first 
recorded by receivers 4 of each well 3. The waves reflected by the discontinuities of the 
zone (seismic interfaces) are propagated upwards. These upgoing waves 10 are also 
recorded by the various receivers 4. The upgoing and downgoing waves are thus 
superimposed in the seismograms. The records are subjected to the conventional 
processings known in the art suited for those obtained by means of the VSP technique 
(Vertical Seismic Profiling). Finally, the respective contributions of the seismic sources to 
the composite vibrational signal are discriminated and the seismograms equivalent to those 
that would be obtained by actuating the seismic sources separately are reconstructed. 
Sinusoidal signals of different firequencies, in their fundamental components as well as in 
their respective harmonics, or signals based on wavelets, Legendre polynomials or random 
series, etc., are for example used as orthogonal signals. 

[0010] In the case notably where the orthogonal signals emitted are sinusoids, 

discrimination of the respective contributions of the seismic sources is for example 
performed by determining the amplitude and the phase of the composite vibrational signal 
at the fiindamental firequencies of the control signals applied to the seismic sources. This 
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discrimination comprises for example weighting of the signals recorded by a bell 
weighting (or apodizing) factor and determination of the amplitude and of the phase of the 
composite signal, or selection by Fourier transform of lines of the complex spectrum 
associated respectively with the various weighted signals. Reconstruction of the 
seismograms corresponding specifically to the various seismic sources is carried out for 
example by applying, after separation thereof, an inverse Fourier transform to the lines 
associated respectively with the various weighted signals. 

[0011] Furthermore, it is well-known that reservoir development induces 

mechanical stress variations in the reservoir and in its overlying bed. These stress 
variations have the effect of generating seismic signals that are propagated in the 

formations. The receivers of the various antennas 2 record these signals (upgoing 
waves 11). 

[0012] Since seismic sources 5 are arranged much closer to receivers 4 of 

antennas 2 than to the reservoir and since the energy that is emitted is much more intense 
than the energy of the induced microseismic signals generated in the reservoir, the 
microseismic signals are generally concealed and undetectable. 

SUMMARY OF THE INVENTION 

[0013] The method according to the invention allows separation of induced 

microseismicity signals and multisource seismic signals within the context of monitoring 
operations in an underground zone under development. 

[0014] It comprises emission of seismic waves in the formation by coupling 

therewith one or more seismic sources. In the latter case, these one or more sources emit 
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simultaneously orthogonal signals so as to form a composite vibrational signal, reception 
of the signals reflected by the formation in response to the emission of seismic waves, 
recording of the signals received by at least one seismic pickup and processing of the 
recorded signals so as to separate the respective contributions of the at least one seismic 
source to the signals received and to reconstruct the seismograms equivalent to those that 
would be obtained by actuating the seismic sources separately. 

[0015] The induced microseismicity signals (passive seismic) are separated from 

the seismic signals emitted essentially by isolating the contribution of each source by 
comparison with a reference spectral model at the emitted fiindamental frequencies and at 
the respective harmonics thereof, and by reconstructing in the time domain the induced 
microseismicity signals and the seismic signals emitted. 

[0016] The spectral contribution of the microseismicity signals to the spectrum of 

the signals received is obtained for example by subtracting the amplitude and phase values 
associated with the reference spectral model from the amplitude and the phase of the 
spectrum associated with the records. 

[0017] According to an implementation mode, a current model formed by updating 

a spectral model taking account of the contribution of previous recording cycles is for 
example used as the reference spectral model. 

[0018] The current spectral model can be formed by determining a mean value of 

the frequency spectra formed from (earlier and/or later) records obtained for the same 
source and the same frequencies. 
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[0019] According to another implementation mode, the current spectral model is 

formed by determining a medism value of the frequency spectra formed from earlier 
records obtained for the same source and the same frequencies. 

[0020] According to another implementation mode, the current spectral model is 

formed by extrapolation or interpolation from the frequency spectrum, from close spectral 
values. 

[0021] According to another implementation mode, the current spectral model is 

formed by using the line obtained at this frequency. 

[0022] According to an implementation mode, discrimination of induced 

microseismicity signals from among signals emitted within the context of active seismic 
monitoring operations is obtained by means of the following stages : 

a) for each recording n of a recording cycle p, respective contributions of the various 
sources at the fundamental frequencies are calculated; 

b) a ratio of a contribution to a current spectral model formed by updating a previous 
spectral model from frequencies emitted during the previous recording and from their 
harmonics is then calculated; 

c) a part of the recording n of cycle p that can be associated with the active seismic 
monitoring operations is deduced; 

d) a part of the recording n of cycle p strictly relevant to the passive microseismic 
activity is deduced; 
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e) seismograms that can be associated with the active seismic monitoring operations 
are formed by inversion in a time domain of respective spectral contributions of each 
seismic source at the fundamental frequencies and at their harmonics thereof; and 

f) the underl3dng microseismic signals contained in the records are formed by 
inversion in the time domain from the part relevant to the passive microseismic activity. 

[0023] The spectral contribution is obtained by multiplying a transfer function 

between a wavelet characteristic of the source and a seismogram associated with a given 
receiver by a wavelet characteristic of the source. This transfer function can be 
continuously updated during a current cycle from an estimation made during a previous 
cycle and from an initial estimation made during the current cycle by the relation. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0024] Other features and advantages of the method according to the invention will 

be clear from reading the description hereafter of an embodiment given by way of non 
limitative example, with reference to the accompanying figures wherein : 

[0025] Fig. 1 diagrammatically shows the monitoring device installed in the field 

for seismic monitoring of an underground formation; 

[0026] Fig. 2 shows a mode of layout of the emission-reception units in the field; 

[0027] Fig. 3 diagrammatically shows one of these emission-reception units 

comprising a seismic source and a reception antenna; 

[0028] Fig. 4 shows a variant of the monitoring device of Fig. 1 ; and 
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[0029] Fig. 5 shows examples of signals and frequency spectra mentioned in the 

description. 

DETAILED DESCRIPTION OF THE INVENTION 

[0030] The method of the invention allows carrying out seismic monitoring 

operations in an underground zone by means of a series of seismic pickups and one or 
more vibrational seismic sources. If several simultaneous sources are used, they are 
actuated simultaneously by signals at different frequencies selected so as to allow 
discrimination of the contributions of each source in the seismograms obtained from the 
signals received and recorded. 

[0031] This is generally performed by control of the various sources by sinusoids 

of different frequencies and by means of known numerical calculation techniques such as 
using inverse Fourier transform to separate the contributions to the seismograms obtained 
from the various seismic sources. 

[0032] Distortion cannot be disregarded with mechanical sources. Simultaneously 

with a frequency fj, a single source Sj emits frequencies 2fi, 35 ... nfj. Consequently, if fi 
and fj are the respective frequencies of the two sources Si and Sj of the series of sources, 
we must have fj fj, as well as ^ 2fj, fi ;t 3^, ... fi nfi in order to be able to separate 
their respective contributions. 

[0033] As already described in the aforementioned French patent -2,805,051, if the 

composite signal Pt consisting of the sum of N sinusoids {fi, Aj, Oi} with 1 < i < N is 
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emitted, all the frequencies being contained in a spectral band contained between two 
boundary frequencies ft and ff, the seismogram T, observed at point R will have as Fourier 
transform at the frequency fi, the number of module Aj and of phase equal to the 
amplitude and to the phase of sinusoid Tj. It is thus possible to reconstruct by inverse 
Fourier transform seismogram by emitting successively all the sinusoids of 
frequencies t to ff. 

[0034] In the case where, for example, all the amplitudes Ai are 1 and all the 

phases are zero, the signal P, obtained is very close to the signal resulting from 
autocorrelation of a sliding-frequency signal contained in the sweep range [fb - fj, 
commonly used in vibroseismic prospecting. According to the discrete Fourier transform 
theory well-known in the art, if one wants to listen to source Si during the time te, the 
frequency increment between the sinusoids is Af = 1/te and the number of sinusoids 
necessary is Nf = (ff - fb)tc. 

[0035] N vibrators SI , S2, Si,. . . Sn installed in the field can thus be simultaneously 

excited by means of vibrational signals whose frequencies are such that each source is 
excited successively by each above sinusoid Nf at any time, on condition that the 
respective frequencies of the sinusoids emitted at the same time by the various seismic 
sources are all different from one another. Nf recordings corresponding to the Nf 
frequencies selected in the sweep range are thus performed with each vibrator. Separation 
of the signals received by the pickups in the field, in response to the simultaneous 
emission of the various signals, is thus obtained by selection of the line at the suitable 
frequency. 
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[0036] Cycles of Nf recordings are thus repeated. The record En of a given 

acquisition cycle contains a set of fundamental and harmonic frequencies which only 
depends on n. Starting from this observation, the method that will be applied to detect the 
effects of microseismicity is based on the observation which has proved true in practice 
according to which the amplitude and phase variations observed at the frequency f and at 
the time t are very little different from those observed at the frequency f - Af and at the 
time t - At, This allows forming a cumulative spectral model constructed from records of 
previous cycles at the same frequencies but formed over a long time during which the 
characteristics of the seismic signals emitted are likely to vary, and permanently adapted 
from records obtained during a sufficiently short period for the variations of the signal to 
be disregarded. The efficiency of method is optimum when the successive frequencies 
emitted by the same source are regularly incremented. Cycles in upgoing and downgoing 
frequencies can for example be alternated. 

[0037] The procedure for showing in practice this spectral contribution of the 

underlying microseismic signals is as follows. 

Notations : 

■ fi,„ designates the fundamental frequency emitted by source i during the recording n 
(this frequency is independent of the cycle). 

In order to simplify the notations of the following quantities, the subscript giving the 
number of the receiver is left out. Whatever this number, the quantities are calculated in 
the same way. 

■ Rp^ designates the recording of rank n of cycle p, 
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■ Cp.i^ the contribution of source i to the recording n of cycle p, 

■ Mp^ the initial model of the active part of the recording n of cycle p. This model 
contains only the frequencies emitted and their harmonics during recording n, 

■ Ep,„ the ratio between the active contributions of the recording n of cycle p and the 
initial model of the active part of the recording n of cycle p, 

■ Ap^ the active part of the recording n of cycle p, 

■ Pp.„ the passive part of the recording n of cycle p (microseismic), 

■ k the differences updating coefficient (typically of the order of 15 to 25 %), 

■ h the models updating coefficient (typically of the order of 5 to 10 %). 

[0038] Considering these notations, the method first consists, for each recording n 

of cycle p, in calculating the respective contributions Cp,,,„ of the various sources at their 
fundamental frequencies and at their harmonics by means of the method described in the 
aforementioned French patent 2,805,051. Owing to the possible presence of microseismic 
energy, the contribution of the various sources is likely to be affected by noise. 

[0039] This estimation can be improved by associating the contributions of all the_ 

receivers of a same recording Cpj^. These contributions can be considered to be the 
product of two factors among which one depends only on the source (and not on the 
receiver) and the other on the receiver, and not on the source. 

[0040] If subscript r represents the receiver, 

Cp,!^^ — Sp,!^ Tp^i^^. 
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T only depends on p, i and n by the variation of the characteristics of the receiver and 
of its environment with time. These variations are very slow and can be disregarded or 
estimated more precisely. 

[0041] S and T are determined to within one factor. The average of the downgoing 

waves can thus be selected for Sp.i^. This average is obtained by means of a calculation 
well-known in VSP processing (flattening of the downgoing waves and average). The 
Tp.i.„^ are then the transfer functions between the downgoing wave and the seismogram. 

[0042] In cases where the variation of T is to be taken into account, updating 

coefficient h can for example be used to update Tp.i^., from Tp,i^.,^ as follows : 

[0043] If Copj^ r designates the initial estimation of contribution C and TOp,,,„,, that 

of the receiver factor, S is calculated by flattening and averaging of the downgoing waves, 
then TOpj.„^ is calculated by the relation : 

Co . 
To . = — ^ 

[0044] The final estimation of the receiver factor is obtained by : 

Tp.i^.r = (1 - h)Tp,i,„.,., + hTOp.i^.r 

and the final contribution of source i to cycle p of recording n by : 

Then the ratio Ep,„ of the current recording to the spectral model by the relation is 
calculated: 
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[0045] This calculation requires knowledge of Mp^ whose calculation described 

below is carried out during cycle p-1. It comprises smoothing the previous frequencies so 
as to limit the effect of the microseismic energy in its estimation. The proposed formula 
only uses the previous records and thus allows real-time conditions. In difficult cases, this 
sum can be replaced by an average or a current median. Deferred-time conditions are then 
necessary (the time required to carry out the second half of the recordings used in 
calculation of the average). 

[0046] If Mp+,^ represents the initial model of the active part of the recording n of 

cycle p+l, then : 

[0047] This model essentially weights an average of the current contribution and of 

the previous contributions. The formula used allows variation of the model : the 
contribution of cycle p-q to the model used for cycle p is multiplied by coefficient (l-h)*» 
which tends to zero all the faster as h is close to 1. This allows model Mp^ to follow the 
seasonal variations, and ratios Ep^ to reflect only the daily variations. Knowing spectral 
model Mp^ and ratio Ep.„ of the contribution to the previously established model, part Ap.„ 
of the recording n of cycle p relevant to the active seismic part is deduced: 

Ap^ = Mp^ Ep^. 

[0048] Similarly, part Pp,„ of the recording n of cycle p strictly relevant to the 

(passive) microseismic activity is calculated by the relation : 
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[0049] The respective spectral contributions Cp.i^ of each source Si at the 

fundamental frequencies and their harmonics, after completion of a cycle of 
measurements, provide, by inverse FFT, the seismograms (in the time domain) used 
within the context of active exploitation. 

[0050] The estimation of the underlying microseismic activity contained in 

record Rp.„ is obtained by inverse FFT from the part Pp.„ relevant to passive contributions. 

[0051] Updating of ratio Ep.„ between the current contribution and the spectral 

model Mp^ is performed in a relatively fast way so as to allow fast variations of the signal. 
Updating of the initial model is performed in a much slower way; in fact, the initial model 
might not be updated. 

[0052] The interest of these updatings is to optimize estimation of the active signal 

when its variations cannot be controlled (seasonal and daily variations). In the most 
favourable cases, they are not necessary and estimation of the active part amounts to the 
initial model. 

[0053] Another possible simplification disregards the seismic noise at the emitted 

frequencies (and, if need be, at the harmonics of these frequencies). It no longer takes 
estimation of the active part and merely cancelling or interpolating the lines corresponding 
to the frequencies emitted (and, if need be, those of their harmonics) is utilized in order to 
obtain the passive part. 

[0054] The harmonic lines of the frequency spectrum of each signal recorded are 

first suppressed, and the amplitude and the phase are interpolated at the line values at the 
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frequency of the source; the record that has not been assigned by the emissions of the 
sources possibly containing a microseismic signal is then reconstructed by inversion. This 
implementation mode is rather suitable in the case of a limited number of sources and/or 
sources having a high signal variability with time. 
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